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ABSTRACT 

It has been suggested that both X-ray bursters and 
millisecond radio pulsars may be strange (quark) 
stars, rather than neutron stars. Confirming (or re- 
jecting) this suggestion may require knowing what 
role strong-field effects of general relativity play in 
the accretion flow of the compact X-ray source in 
low-mass X-ray binaries (LMXBs). We discuss the 
range of rotational and orbital frequencies, and of 
masses expected in various models of strange stars, 
and compare them with observational constraints, 
suggested by the observed frequencies of kHz QPOs. 
We explain why future observations of transients 
(such as SAX J 1808.4-3658) may be crucial to under- 
standing the precise nature of the accreting source. 
For flattened (e.g., rapidly rotating) distributions of 
matter, an innermost (marginally) stable orbit may 
be present even if relativistic effects are negligible. 
Depending on the stellar rotation rate, the same 
value of orbital frequency in the innermost stable 
orbit (say, 0.9 kHz) can correspond to a star of mass 
equal to 2.4M , 1.4M , 0.1M©, 0.01M Q , or less. 

Key words: dense matter- equation of state - stars: 
binaries: general -X-rays: stars. 



1. INTRODUCTION 

Bodmer (1971) suggested that quark matter (com- 
posed of roughly equal number of up, down, and 
strange quarks) is the ground state of matter. 
Brecher and Caporaso (1976), Witten (1984), Alcock 
et al. (1986), Haensel et al. (1986), and others, have 
computed the structure of gravitating objects com- 
posed of such matter, the so called strange stars. 
Because the density of quark matter would exceed 
only slightly nuclear density, the radii of these quark 
stars would be comparable to the radii of neutron 
stars, as would be the maximum masses — it would be 
hard to distinguish the two at first sight. Although 
young glitching radiopulsars are likely to be neutron 
stars (Alpar 1987), we are unaware of any reason 
why either the observed millisecond pulsars, or X- 
ray bursters (as well as the "Z sources") in LMXBs, 
should not be quark stars instead. 



Unlike for neutron stars, which can exist only above 
a certain mass (~ 0.1M©), there is no lower limit 
to the mass of a strange star, so the discovery of 
a star with a minute radius, less than about 5 km, 
or of a heavy planet of an enormous density, would 
immediately confirm the existence of quark matter. 
An unusually high rotation rate could also reveal a 
quark star. 

A precise determination of the mass is in princi- 
ple possible through the measurement of the or- 
bital period in an accretion disk (as in the super- 
massive black hole candidate NGC4258, where the 
velocity of water masers at various distances from 
the central object has been measured: Miyoshi et 
al., 1995). In LMXBs, it is currently not feasi- 
ble to resolve spatially any "blobs" orbiting in the 
disk; however, the mass of the compact star can 
be determined from the maximum orbital frequency, 
forb,max, on the assumption that this maximum is 
attained in the marginally stable orbit of general rel- 
ativity (Kluzniak et al. 1990). The recently discov- 
ered kHz QPOs can in principle be used to determine 
the orbital frequencies. In the simplest model the 
QPO frequencies are the orbital frequencies and, on 
the assumption just stated, the mass of the compact 
star in some sources can be thus determined (e.g. 
Kaaret et al. 1997, Zhang et al. 1998, Kluzniak 
1998, Bulik et al. 1999b). Such considerations make 
important determining the range of possible orbital 
periods for various models of quark stars (Stergioulas 
et al. 1999; Zdunik et al 2000a, b) and of neutron 
stars (Thampan et al. 1999). 

Here, we summarize our studies of quark star masses 
and frequencies (rotational and orbital), and point 
out a new complication: the marginally stable or- 
bit exists also in Newtonian gravity! If higher mass 
multipolcs are present in the distribution of matter, 
estimates of the central mass from the maximum or- 
bital frequency alone are unreliable. As we show, the 
same value of forb,max can correspond to two stars 
differing in mass by a factor of a hundred, or more. 
It is crucial to know the rotational period of a quark 
star, if its mass is to be even approximately deter- 
mined from the orbital frequency. 
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Figure 1. Maximum orbital frequency for test par- 
ticles in circular motion in the equatorial plane of 
a quark star with baryonic mass Mb — 0.01M Q; 
as a function of stellar spin rate (thin continuous 
line). At low rotation rates, Keplerian orbits graze 
the stellar surface, but for f rot > 970 Hz a gap sep- 
arates stable circular orbits from the stellar surface, 
and the maximum orbital frequency is attained in the 
marginally stable orbit (ISCO). Also shown are the 
ISCO frequencies for maximally rotating quark stars 
(thick continuous line, from Stergioulas et al. 1999). 
All calculations were performed with fully relativis- 
tic codes for the MIT-bag model, with a = 1/3, po = 
4.28 x 10 14 g/cm 3 ; in eq. [1]. The dotted lines are 
extrapolations of these calculations to lower orbital 
frequencies. For maximally rotating models, f rot is a 
monotonically increasing function of mass, all stars 
on the thick continous line have M > 1.3M . Thus, 
the same value of f or b,max (e.g., 1.1 kHz) can corre- 
spond to the ISCO around a 0.01M Q star rotating at 
1kHz, or a 2M Q star rotating at about 1.3 kHz. 



2. TWO MODELS OF QUARK MATTER AND 
OF STRANGE STARS 

Traditionally, strange stars were modeled with an 
equation of state based on the MIT-bag model of 
quark matter (Farhi and Jaffe 1984), in which quark 
confinement is described by an energy term propor- 
tional to the volume. Recently, Dey et al. (1998) 
proposed a model in which quark interactions and 
confinement are described in a more self-consistent 
manner. In fact, both models give a simple equation 
of state 

P = a(p- Po )c 2 , (1) 

where a is a constant of model-dependent value (close 
to, but generally not equal to, 1/3 for the MIT-bag 
model, Zdunik 2000). It is the value, p , of density 
at zero pressure which is crucial for various limiting 
properties of strange stars; e.g., the maximum mass 
of static strange stars is M = 2M Q (p /A x 10 14 )~ 1/2 
in the MIT-bag model (Witten 1984), and the scaling 

of mass with Pq 1 ^ 2 is general. 

Below, we describe our main results obtained for 
MIT-bag model stars. For other models of quark 
stars the results are qualitatively similar, However, 
in general the Dey et al. (1998) model gives more 
compact stars, with lower maximum mass and lower 
radii than the MIT-bag model stars. This also trans- 
lates to higher possible rotational and orbital fre- 
quencies for the "Dey stars" (Gondck-Rosihska et 
al. 2000a,b). 



3. SUMMARY OF RESULTS 



We were mostly interested in the value of maximum 
orbital frequency of test particles (or fluid elements 
in an accretion disk) around such stars. A very 
specific question was asked (Bulik et al. 1999a, b) 
whether the observed maximum frequency of kHz 
QPOs in 4U 1820-30, 1.07 kHz, is compatible with 
strange star models — for static stars such a frequency 
is obtained in the marginally stable orbit only for 
the relatively high mass of nearly 2M©, according to 
the formula for the maximum orbital frequency in 
Schwarzschild metric, f ms = 2198kHz(M /M). We 
found, that with the current physical constraints on 
po (eq. [1]) masses of static strange stars as high as 
2.4M Q cannot be ruled out (Zdunik et al. 2000a). 

To gauge the range of allowed orbital frequency, Ster- 
gioulas et al. (1999) computed models of maximally 
rotating quark stars. The main results are, that for 
these models the innermost (marginally) stable cir- 
cular orbit (ISCO) is always well above the stellar 
surface (even though rapid rotation leads to a large 
increase of the equatorial radius), and that the ISCO 
frequencies are usually much lower for rapidly rotat- 
ing stars than for static ones, with a value as low 
as 950 Hz for a 1.4M Q star rotating at the equato- 
rial mass shedding limit. Only close to the maxi- 
mum mass do ISCO frequencies in maximally rotat- 
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ing quark stars become comparable to the static val- 
ues. However, the ratio T/W in these stars exceeds 
the Newtonian limit for instability (T/W > 0.1375), 
so probably these extreme configurations cannot ex- 
ist. Detailed models of quark stars rotating at all 
possible rates can be found in Zdunik et al. 2000b. 

The gap between the ISCO and the stellar surface 
is always present for quark stars more massive than 
1.4 to 1.6M Q (depending on the model). At some- 
what lower masses the gap disappears for moder- 
ate rotation rates (but the marginally stable orbit 
is present for both slowly and rapidly rotating stars 
of this mass). It is only for the lowest mass stars, 
that the gap disappears even at zero rotation rate. 



4. THE LOW MASS LIMIT 



The limit M — » corresponds to p — > p in cq. (1), 
i.e., the star tends to a nearly uniform density con- 
figuration. At the same time, the dimcnsionlcss an- 
gular momentum parameter j = cJ/(GM ) — > oo, 
so an expansion in powers of angular momentum 
(Kluzniak and Wagoner 1985, Shibata and Sasaki 
1998) is impossible; however, the limit is, in fact, 
Newtonian. In Fig. 1 we show the maximum value 
of orbital frequency, as a function of stellar rotation 
rate, for a low mass quark star (the curve is nearly 
universal for all M < 0.1M Q ). The computation was 
carried out with a fully relativistic code, described in 
Gourgoulhon et al. (1999). 

For all periods of rotation larger than about 1.03 
ms, there are stable circular orbits just outside the 
surface. For non-rotating stars, f rb,max is simply 
given by Kepler's law, and its value decreases as the 
equatorial radius increases for rotating stars. The 
cusp in the curve at about 970 Hz corresponds to 
the appearance of a gap between the star and the 
ISCO. In this case, this is a purely Newtonian effect 
related to the flattening of the star at rapid rotation 
rates (see also Zdunik and Gourgoulhon 2000). 

To demonstrate the existence of marginally stable 
orbits in Newtonian theory, we present in Fig. 2 
the angular momentum in circular orbit around a 
uniform density disk, computed from the Newtonian 
gravitational potential of this non-spherical mass dis- 
tribution. 

The conclusion is, that because of the appearance 
of the marginally stable orbit in Newtonian gravity 
(when octupole and higher mass multipoles are suf- 
ficiently large), the maximum orbital frequency can 
reach fairly low values even for extremely low-mass 
stars. Conversely, / = 1.1 kHz in the marginally 
stable orbit need not imply that M = 2M Q , we 
have exhibited a specific example where the mass is 
< 0.01M©, instead. 

This work has been funded by the KBN grants 
2P03D01816, 2P03D00418. 




Figure 2. Specific angular momentum squared, L 2 , in 
Newtonian theory of circular orbits around a homo- 
geneous disk of radius R and mass M, as a function 
of the orbital radius r (continuous line). Also shown 
are succesive approximations to this angular momen- 
tum, when the multipole expansion of the Newtonian 
potential is truncated at the monopole term (dotted 
line), and at the quadrupole, octupole, 16-pole terms 
(successively higher curves). Note that for this disk, 
in purely Newtonian gravity, the marginally stable 
orbit — which corresponds to the minimum of L 2 — 
appears (just outside the disk) already in the octupole 
approximation, but actually occurs at r w 1.3i2, with 
the corresponding orbital frequency, 2tt forb.max — 
L/r 2 , about 20% higher than for a circular orbit 
(of the same radius) around a spherical mass dis- 
tribution of mass M. On the other hand, the same 
forb,max is only 0.8 of the value of the orbital fre- 
quency at r = R around a spherical mass M. 
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